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1 Abstract 

A Dynamic Programming based polynomial worst case time and space algorithm is described for computing Hamil- 
tonian Path of a directed graph. Complexity constructive proofs along with a tested CH — h implementation are 
provided as well. The result is obtained via the use of original colored hypergraph structures in order to maintain 
and update the necessary DP states. 

2 Introduction 

Dynamic Programming hardly requires introduction. Since the term was in introduced by Richard Bellman in 1940's 
it has countless applications and it's power allows to have fast algorithms in cases where at a blush there seems to be 
no way to avoid exponential time. Such examples include graph algorithms: Bellman-Ford, Floyd- Warshall as well 
as option pricing, numerical solutions to HJB equations ("backwards in time"), discrete optimal control policies, 
knapsack small block size problem [2], Smith- Waterman local sequence alignment algorithm and many others. 

Hamiltonian path (H-path) in directed and undirected graph is one of Karp's 21 NP-complete problems [3]. 
The question it asks is to find cycle or path in a given graph which visits every vertex exactly once. Therefore one 
can see that H-path is the shortest path which visits all nodes in the graph and thus provides the answer to the 
transportation problem in question. 

David Zuckerman [4] showed in 1996 that every one of these 21 problems has a constrained optimization version 
that is impossible to approximate within any constant factor unless P = NP, by showing that Karp's approach 
to reduction generalizes to a specific type of approximability reduction. Significant progress has been made in 
expanding the class of graphs for which polynomial solution does exist: notably Ashay Dharwadker [S] discovered 
an algorithm for a broad class of highly connected graphs. 

The intuition behind the presented approach comes from physical objects that seem to have a way to solve 
some NP-complete problems such as soap bubbles forming (almost) minimal surfaces and protein folding satisfying 
hydrophilic-hydrophobic Boolean conditions. 

Consider a cobweb attached to a tree with one end and being pulled by the other end. One can observe that the 
lowest (average) strain or highest slack in the cobweb is achieved along the longest sequence of web segments. This is 
of course a Hamiltonian path in the graph represented by the cobweb. If one wants to find it then it makes sense to 
implement some incremental scheme which akin to a difference scheme for PDF with boundary condition, computes 
the strains incrementally starting with one attachment point and moving towards the opposite attachment point. 

3 Definitions 

Denote G,V{G), E{G) - the input directed graph, it's vertices (aka nodes) and edges respectively, n — \V{G)\. 
Assume G to be a general directed graph which does not have multiple edges between any two vertices and that 
there are no loop edges. 

A given instance of H-path problem is represented by a triple (G, s, e) where s, e € V{G) and the task is to find 
path 7r(s, e) G V{G) starting at node s and ending at e so that every vertex in V{G) is visited and only once. 

In the following description I shall opt for a less conventional but more compact and practical terminology in 
describing the main algorithm by adopting Object-Oriented terminology. The entities used here will be Objects 
which consist of Attributes - representing data and Methods - representing functions as well as traditional program 
functions and variables - these are distinct from object methods and attributes in that they are "static", effectively 
having global scope and accessible from all objects. 
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If an object A is an attribute of _B I typically write As unless making a statement about any object A or it is 
clear from the context what B is. Also instead of single letter notation I elect mnemonics that are much shorter 
than full object names but are easier to associate with full object names despite sounding weird. 

I start with introducing a concurrent graph traversal method which unlike DFS/BFS (Depth First Search/Breadth 
First Search) has some useful properties in this context. This will become a basis for the workflow of main algorithm. 

Denote P{G) family of H-paths in G between nodes s, e. 

Definition 3.1. Let botmaster (C++ code: edg) be an object with the following: 

• Attributes: 

(1) bots - family of all bot objects, defined below. 

(2) march step marchStep G l..n 

(3) dock a Boolean variable associated with pair (v , marchStep) , v € V{G) 

• Methods: 

(1) march for all bot G bots invokes advanccbot, defined below. 

(2) addDelete{vi,V2) (C++: edg::add_delete) defined below 

Definition 3.2. Let bot be an object invoked by master object botmaster, which on march step marchStep 
conditionally transitions from node a to its neighbor b along edge ab, denote that fact as bot{a, s) — >■ b. Bot's only 
attribute is current node a e V{G) and it has 

• Methods: 

i advance: 

(1) bot{a, s) ^ b ^ dock{b, s) — F, dock(a, s — I) = T. This means dock acts like mutex for bots - first 
bot accessing b docks at b. 

(2) addDelete{a, b)},otmaster is invoked 

(3) If bot{a, s) — > bi...bh, h > I, bot{a, s) spawns botlets bot{bi, s + l)...bot{bh, s + 1) 

ii terminate - deletes bot object under the following conditions: 

(1) If ^b : bot{a, s) — > 6, bot(a, s) is terminated. 

(2) bot{a, n — 2) ^ e and there are no other v : bot(a, n — 2) v hot march stops. 

Definition 3.3. Let hist{bot) be the history of nodes visited by bot, this includes "genetic memory", i.e. newborn 
botlet inherits history of it's parent. 

Remark 3.1. There are no more than n bots at any given time. 



This follows from correspondence between bots and nodes guaranteed by dock condition (3.2). 
Remark 3.2. Let tt is H-path on G. Then there exists a bot: hist{bot) ~ n. 



Proof. Assume that tti — (s, ei) C vr, is the longest path traversed by any bot in |7ri | steps. By (3.2 1.1 we have a 
bot docked at ei which will traverse to 62 : (ei, 62) C tt on step |7ri| + 1 thus violating the assumption of maximality 
ofTTi. □ 

Remark 3.3. Any vertex of G is visited at most n times. 
This follows from condition (3.2).ii.(2) 

Definition 3.4. Let 

• Color G of graph g for node w G G be an integer corresponding to a class of paths {(s...w)} in graph g. 

• colors{g) be all colors of graph g. Also: 

• Each color has a unique node it relates to via surjective function: cono : conoiC) — >■ V{g) (stands for color 
nodes: C++ dg::color_nodes). 

• Denote 7r(C) set class of paths for color G. Note that |7r(C)| may be exponentially big and therefore this is 
not explicitly computed by the algorithm. 
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• There is a base color which s is painted in - emptycolor 
Definition 3.5. Let, 

• cohi (color hierarchy, C++ color_hierarchy) be a directed graph of colors, where edge ciC2 denotes that every 
path over nodes of g in ci is a subpath of some g-path in C2, has form {(c, Sc), Sc C colors{g)} 



colors{cohi) = V{cohi) returns colors of cohi 



• sucoiC) is a global function (C++ sub_colors) which returns all colors in cohi which have path to C, i.e 
sub-colors of C , also has form 



{{c,Sc),ScC. y color s{gpagra(v, slack))} 

veV(G),slack=l..n-l 



• noco(v) (node colors, C++: dg::node_colors ) is colors over node v e V{g)^ has form {(v, 8^)^ C colors{g)} 

• cn{v) (colors for node, C++: dg::cn ) is all colors C over v such that cono{C) — v 

• top{cohi) , base{cohi) are respectively first (no ancestors) and last (no descendants) colors in cohi graph. 
Lemma 3.1. Vu G g, \cn{v)\ = 0{n) 



Proof. Be remark (3.3) we have 0{n) bit visits to v. Each visit performs merge for given pagra{v, slack) which 
generates new C : cono{C) ~ v. □ 

Definition 3.6. Color C is called inactive: inactive{C) iff 

• cono{C) = 

• or cono{C) ^ V{g) 

• or noco{cono{C)) — 

• or C ^ top{cohi) A cohi{C) — 

Definition 3.7. Let S £ colors (cohi), a e V{g) dep{S\cn{a)) is a subset of S such that either 

C g 5, ^Tr{C,top{cohi)) £ cohi : cn{a) C] Tr{C,top{cohi)) — 

Or 

C £ S, ^n{C, base{cohi)) £ cohi : cn(a) H 7r(C, base{cohi)) = 

Definition 3.8. For all a £ V{G), slack £ l...n call path graph pa(7ra(a, slack) (C++: dg) an object which includes 
all paths (s...a) C G of length slack and has: 

• Attributes: 

(1) Directed graph g with nodes in V{g) 

(2) Node colors noco{V{g)) = {{uk -> {c^, .., cjj})}, noco is a one-to-many map V C 

(3) Color nodes cono{C) 

(4) Color hierarchy cohi{g) 

• Methods: 

(1) addSlack(b), b ^ V{pagra{a, slack)), does: 

i adds edge ah to gs , Vm 

ii invokes paint which performs: 

(1) allocates new C 

(2) Vw £ V{g),noco{v) := {noco{v),C) 

(3) cono{C) := b 

(4) cohi{C) :— top{cohi) 
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(4) suco{C) := colors{cohi) 

(2) Remove node (C++: dg::rni_node) reno{a): 

i deletes noco{a) 

ii executes bleach{cn[a]) which: 

(1) computes D = {dep{cn[v\\cn[a])\iv € V{g)} 

(2) VC e D remove C from all attributes of pagra{a, slack) 

(3) ensures VC e color s{cohi) , inacfAve{C) = F 

(4) every color C € co/ii has a direct predecessor or bleach removes such C 

(5) if noco{v) — 0, remove v from g \/v ^ ^io) 

(3) Merge pagra{a) + pagra{b): unions are taken for cohi,cono,noco, i.e. it has form: 

{(a, Sa)} + {{b, Sb)} = {a,b: a = b^ {a,SaU Sb), otherwise{a, Sa), {b, Sb)} 

Definition 3.9. Denote 

sucoDFS{a,X,suco,dir),a G color s{cohi),X c color s{cohi),dir G {up, down} 

a.pagra method (C++: dg::has_path_to_top()) which traverses cohi starting from color a, bypassing colors X (C++: 
dg::xcolors) and uses suco to make transitions between nodes as follows: 

(1) dir = up 

(1) run regular DFS against graph cohi starting at color a 

(2) when next node is c G X return to the next available sibling node 

(3) transition from color vertex ai to 02 only if tti C suco{a2) where tti = (a...ai) is a sequence of colors - 
path - in cohi traversed before 02- 

(4) terminate and return path n when next node is in top{cohi) 

(5) when all colors have been traversed terminate and return 

(2) dir = down 

(1) run regular DFS against graph inv{cohi) starting at color a 

(2) when next node is c G X return to the next available sibling node 

(3) transition from color vertex ai to 02 only if Vc € 7ria2 G suco{v) 

(4) terminate and return path tt when next node is in base{cohi) 

(5) when all colors have been traversed terminate and return 

Remark 3.4. By definition of sucoDFS we have: 

sucoDFS{a, X, suco, up) = V sucoDFS{a, X, suco, down) = dep{a\X) = {a} 
Definition 3.10. Given directed graph G inv{G) inverts arc orientation. 
Definition 3.11. Slacks slacks{a) is an object, effectively a wrapper for pagra{a, slack) with: 

• Attributes: 

(1) A single node v gV{G) 

(2) A set of slacks S 

(3) Path graphs pagra{a, slack), slack G S 

• Methods 

(1) Add slack - addSlack{b),b ^ V{gpagra(a,siack)), invokes addSlack{b)pagra(b,siack),'^slack 

(2) Remove node: reno{a) invokes reno{a)pagra(b, slack) ^^slack 

(3) Merge slacks{a) + slacks{b): invokes pagra{a, slack) + pagra{b, slack), M slack G S 
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Definition 3.12. addDelete{a,b)botmaster is a method which performs the following steps: 

• renOsiacks{a){b) 

• addSlacksiacks(a){b) 

• Merge slacks: slacks{b) := slacks{a) + slacks{b) 



Figure 1: Main objects 




Definition 3.13. For a connected sequence of nodes tt dup{n) is a set of repeated nodes. 
Definition 3.14. Denote P'{G) class sucoDFS'-traversible paths Vs, e e V{cohi),X = 0. 

4 Explanation 

The idea here is to categorize families of paths without explicitly describing each one of them. Slacks object for a 
given node encodes paths to a given node over G and groups them by length (slack) . For each given slack pagra 
object is the encoding of paths given by all paths traversible in it's own copy of graph g. Graph g by itself is not 
enough to keep track of those paths due to the way it has to be updated: via merge and reno operation, reno 
operation ensures paths don't run over the same node twice since addDelete removes each node v from g when 
hot{u, march Step) v. After that, merge executes pagra{u, slack) + pagra{v , slack) (ii pagra{v , slack) / 0). If 
91 = 9pagra{u, slack) and ^2 = gpagra{v , slack) both Contain somc nodc w then the resulting g will lose track of paths 
through w resulting in "synthetic path" which is not a path in G for a given slack. To preserve actual paths, concept 
of color is introduced. However just the color alone is not enough since we have cases when we start with some 
color, for example, "red" in slacko which then "splits" into slack\ and slack2 each one carrying it's own version 
of "red" which then gets modified by their own reno operations. After a while slacks combine again but now red 
two distinct path classes. We avoid exponential explosion in number of colors by using cohi, which associates each 
copy of red with some own unique color added in slacksi, for example, "blue" and in slacks2 - "brown". These 
new colors are added as part of addDelete method in bot march and so their number is "small" as shown below. 
When slacki + slack2 is computed, suco includes {brown, {red, ...)), {blue, {red, ...)). Method sudoDFS uses suco 
to avoid "synthetic path" as shown below. 
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5 Path existence 

Lemma 5.1. Total number of colors is \colors{suco)\ — 0{n^). 
Proof. Indeed, a color is assigned by paint method which is called for: 
• every vertex in w € V{G) - 0{n) 



every bot visit to w - 0{n) by Remark (3.3) 
every pagra{v , slack) - 0{n) 



Lemma 5.2. Let tt E P'{cohii + cohi2) then 

/Stti C P{cohii), Tr2 C P' {cohi2) : tt — (7ri,7r2) 
i.e. sucoDFS does not return n that starts in cohii and continue in cohi2. 
Proof. Assume such path tt existed, then we would have 

ci e color s{cohii)\color s{cohi2), C2 £ color s{cohi2)\colors{cohii) : Ci,C2 € n 



□ 



See Figure 2 for illustration. But then we must have ci £ suco{c2) by rules (3. 7). 1-2 which is impossible because ci 
is not even in colors{cohi2) 



Figure 2: Impossible synthetic path on merger 




o 
o 



cohi 1 & 2 
cohi 2 
cohi 1 



□ 

Lemma 5.3. Let j be a sucoDFS-traversible path, tt = cono{-^), u £ t: and importantly v £ g but v ^ n, then 

dep{cn{u)\cn{v)) = 0,Vu G tt 

Proof. By definition of dep, Vit £ tt, cn{u) has path to top{cohi) and that path is 7 and 7 n cn{v) = therefore 
dep{cn{u)\cn{v)) =0. □ 
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Lemma 5.4. Let 

TT = cono{sucoDF S{top{cohipg^gra(b,j)),^, suco, down)) 

then dup{T:) = 0, |7r| = j . 

Proof. Assume there was an intersection dup^ir) ^ then we have ci,C2 : cono{c\) — cono(c2) however if ci 



was traversed after c\ then by (3.7 1.1 observe that c\ G suco(c2) but by suco construction this is impossible if 
cono(c\) = cono{c2). 

\n\ = j folows from bleach method definition: every color C £ cohi has a direct predecessor or bleach removes 
C, also paint method and cono surjectivity ensures that any color path in cohipagra{b,j) from base to top is j long. 

□ 



Lemma (5.4) allows us to recover final H-path from cohipagra(e.n-i) using standard sucoDFS function. That is 



if such path is included in cohipagra{e.n-i)i lot's show this is indeed the case: 
Theorem 5.5. If P{G) ^ then Btt e P{G): 

TT = cono{sucoD F S {top{cohipngra(e,n-i)) J 0; SUCO, down)) 

Proof. Using induction by march step, this trivially holds at node s, marchStep — 1, assume at marchStep = k 
we have: 

7j = sucoDFS{top{cohipagra{vj,j),^T suco, down) 
Let TTj = cono{'jj). When bot(vj,j) — >■ Wj+i we can show that, for ttj+i = {TTj,Vjj^i) 

TTj+i = cono{sucoDFS{top{cohip„_gra{vj+i.j+i)j^T suco, down)) 

Indeed bot(vj,j) — > Wj+i invokes addDelete which invokes methods merge and reno, merge does not affect tt^ 



by lemma (5.2) while if Wj+i G y{gpagra(vj,j)) reno will invoke bleach{dep{cn{v)\cn{vj^i))),yv G tt^ but since tTj is 



sucoDFS traversable and w^+i ^ Hj therefore by lemma ( |5.3[ ) 

color si^j) <f. dep{cn{v)\cn(vjj^i),\lv G tTj 
therefore 7j and therefore tTj is not affected by bot{vj,j) — > Wj+i. □ 

6 Algorithm 

• Start with one bot which current node set to s 

• slacks object and the associated pagra objects are initialized 

— slacks consists of just one pagra 
^ dpagra cousists of Only One node s 

— cohi consists of just empty color 

— noco is just (s, emptycolor) pair 

— cono is (emptycolor, s) 

• Bot advances s — > wi according to it's rules, and invokes addDelete function with updates slacks{vi) 

• Steps above are repeated, until pagra{e,n — 1) G slacks{e), i.e. we have a path graph with slack equal to 
\V{G)\-l. 

• Extract the H-path by computing 

cono{sucoDFS{top{cohipagra(e,n-i)),^, SUCO, down)) 
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7 Complexity 

Time complexity is based on that: 

• there are 0{n) G nodes 

• visited by bots at most 0{n) times 

• each visit requires runo call for each of 0{n) graphs pagra{v, slack) 



reno, requires sucoDFS call for each color a of at most 0{n) sized (lemma (3.1 1) cn{v) color set 



sucoDFS{a,cn{b),suco,dir) complexity is 0{n ) since: 



any sucoZ? i^S'-traversible path in cohi: ■n{hase{cohi)^top{cohi)) is 0{'n) long by lemma (5.4) 

set X ~ cn{b) cardinality is at most 0{n) by lemma ( |3.1[ ) and therefore DFS requires at most 0{n) 
rollbacks when X color is encountered 



— validating conditions ( |3.7[ ).1,2 on every step costs 0{n) 

• which is executed for each 0{n) nodes w e .g to decide if the node is bleached 
Space complexity is based on that: 

• there are 0{n) G nodes 

• each node has 0{n) pagra objects 

• each pagra object has cohi as biggest attribute 

• each cohi attribute has at most 0{n^) colors 

Thus total current worst case time complexity is O(n^) and space complexity O(n^). Average case complexity 
is typically much lower because in practice cn{v) is 0(1) sized and average number of node bot revists is also 0(1). 

sucoDFS{a,cn(b),suco,dir) complexity can be lowered by 0{n) if we eliminate cn{b) from cohi first which is 
done once for all colors a. Than for each a we only have to run sucoDFS{a,(l), suco,dir) in 0{n). 

Complexity can likely be further lowered by an additional 0{n) by consolidating all cohi per G node and across 
pagra of different slack levels. Slack values can be extracted from cohi directly when necessary. 

Additionally dep{cn{b)\cn{a)) can possibly be found faster if we start with an encoding of all swcoDFS'-traversible 
paths P to top. Then on every update of cohi instead or rerunning sucoDFS we could just update P. 



8 Implementation and testing 

The algorithm was tested on 10000 randomly generated graphs each with 17 nodes and 3 outbound edges per node. 
Test graphs were generated using function dg :: gen^graph which is part of the source code and is as follows; 

• Generate a random H-path first. Sample without replacement from a discrete uniform random distribution 
over l...n : U{l...n). 

• For each node on selected H-path generate given number of edges 5 by choosing vertex 1 sequentially from 
the path and then choosing vertex 2 by sampling without replacement from U{l...n) S times. 

Please note that performance tuning was not a priority therefore the author is aware of a number of inefficiencies. 
The algorithm is available for download at 

https://docs.google.eom/folder/d/0B3s6PXhKJO6HWnE4c0VZbVJodDQ/edit . 

It is free to use for academic and educational purposes use but requires a license for commercial and government 
use. Potential commercial and government users should note the algorithm has a Patent Pending status with 
USPTO. 
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Appendix 



Figure 3: Source Code 
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} 

return t n ; 



} 

i nt 



p a i 


nt() { 










i nt 


r et =0; 










per 


m_ p a i n t ( n 


d e_ c 


1 


r s , 


n e w_ c i r ) ; 


coi 


r _ ti i e r a r 


c ti y [ n 


e w_ 


coi 


or] =top_coi ors; 


coi 


r _ n d e s [ 


n e w_ c 


1 


r] = 


:top_node; 


ret 


= ne w_ c i 


r ; 








coi 


or s i a c k s 


[new 


c 


or] 


=5 i a c k i e V e i ; 



} 

voi 



top_coiors.clear(); 

top_coi ors[new_col or]=true; 

subcoi ors[ new_col or] =fl atten(coi or_lii erarchy); 
n e w_ c i r + + ; 
return ret; 



d ciear_inacti 


ve_ 


col 


or s ( ) 


i i St <i nt > r ml 








map2: : iterato 


r i 


t ; 




for ( i t =col or 


hi 


er a 


rchy. 


int c = ( 


F nt 


) i t 


- >f 1 r 


if ( ! acti 


ve 


col 


or(c) 


ma p 1 : : i t e 


rat 


or 


it; 


ma p 1 X = 


col 


°r _ 


hi er a 


for ( j t =x 


. be 


gi n 


0; i 



;gin();it!=color_hierarchy.end();it+ + ) { 
rml.push_back(c); 



=x. end( ) ; j t ++) { 
int sc =(int)jt->first; 
if ( ! a ct i ve_ col or ( sc ) ) { 

col or_hi erarchyl c] . erase! sc) ; 



<i nt >: 


: i 


t e r a 1 r lit; 




( 1 1 t =r 


■ ml 


.begin();lit!=rml.end(); 


1 1 t + + ) 


1 nt c 




( i nt ) * 1 i t ; 




col or_ 


hi 


erarchy. erase(c); 





for 



} 
} 

dg& operator +=(dg& g2) { 
add_cnt + + ; 

merge_maps(col or_hi erarchy, g2. col or_hi erarchy) ; 
merge_maps(col or_nodes, g2. col or_nodes) ; 
merge_maps(descendants, g2. descendants) ; 
j oi n_maps(&top_col ors,&g2.top_col ors); 
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upd_precendants(); 

merge_maps(node_col ors, g2. node_col ors) ; 
get_nodes( ) ; 

clear_inactive_ colors!); 
return *thi5; 

} 

void u pd_ p r e c e n d a nt s ( ) { // from descendants 
map2::iterator it; 
ma pi cd_desc; 
precendants.clearl); 

for ( i t =d e s c e n d a nt s . beg i n ( ) ; i t ! =des c e n d a n t s . end ( ) ; i t ++) { 
int node = it->first; 
mapl &descs = i t ■ >second; 
mapl::iterator jt; 

for ( j t =des c s . beg 1 n ( ) ; j t I =d es c s . e n d ( ) ; j t + + ) { 
int node2=jt->first; 
pr ecendants[ node2] [ node] =t rue; 

} 

} 

} 

ma p 2 1 c h ; 
voi d 



) ; i t I =c i r _ h i e r a r ch y . e n d ( ) ; i t + + ) { 
=0) i c h[ s up_ c i or ] =empt y ma p; 

=descs. end( ) ; j t ++) { 

;t ; 

or ] =t r ue ; 



i nverse_ch( ) { // f r 


om c 1 


r _ 


i ch. cl ear ( ) ; 






map2::iterator it; 






ma pi cd_desc; 






ma p 1 e mp t y ma p ; 






for (it=color_hierarc 


h y . beg 


i n 


int sup_color = i 


t ■ >f i r 


St 


if (ich.count(sup 


coi or 


) = 


ma p 1 Sid e s c s = i t ■ 


>secon 


d; 


ma p 1 : : i t e r a 1 r j t 






for ( j t =des c s . beg 


i n ( ) ; 


i t 


i nt sub_col or 


=j t - >f 


i r 


1 ch[ sub_col or 


1 [ 5Up_ 


c 



} 

ma p 2 c n ; 








void inverse 


_col 


r _ nodes ( ) { 




ma p <i n t , 


1 nt >: 


:lterator it; 




for (1 t = 


color 


_nodes. begi n( ) ; i t 


I =coi or_ nodes 


1 nt 


color 


= ( i nt ) 1 t - >f i r s t 




if ( 


! act 1 


ve_col or ( col or ) ) 


c n t i n u e ; 


c n [ c 


ol or _ 


nodeslcolorll[coi 


or ] =t r ue ; 



end();lt+ + ) { 



} 

} 

void per mpai nt ( map2 &node_ c 1 r s , i n t coior){ 
get_nodes( ) ; 
ma p 1 : : i t e r a 1 r mi t ; 
int node=top_node; 
wh i i e (true) { 

if (node = = s tart_node) break; 

if ( p r e c e n d a n t s [ n d e ] . s i z e ( ) >1 ) break; 

int i = precendants[node].size(); 

node =(int)precendants[nodel. begi n()->flrst; 

} 

for ( mi t =nodes . begi n( ) ; mi t I =no d es . en d ( ) ; ml t ++) { 
int node =(int)(* mi t). first; 
node c ol or s [ node] [ col or 1 =t r ue; 

} 

} 

void get _t op_ c ol or s ( I nt c.mapl &ret) { 
map<int,lnt>::lterator It; 

for ( I t =col or_nodes. begi n( ) ; I t I =coi or_node5. end( ) ; I t + + ) { 
int color = (int)lt->first; 
int node = (Int)color nodes ( c i r ] ; 



APPENDIX 



14 



C:\mydoc5\deform_\sd9.h 



i f 
} 



; node==t op_node] { 
r et [ col or I =t r ue; 



bool active_color(int c] { 

f ( c 1 r _ n d e s . c u n t ( c ) = = 0) return false; 
nt cn=color_nodes[c]; 
f ( n d e s . c u n t ( c n ) = = ) return false; 
f ( n d e_ c 1 r s . c u nt ( c n ) ==0 ) return false; 
f ( n d e_ c 1 r s [ c n 1 . c u nt ( c ) ==0 ) return false; 
return true; 



ma p 1 a I I _ b I e a 

void rm_node( 
g r _ r m_ n d 
ma p 1 : : I t e 
ma p 1 e mp t 
I nverse_c 
ma p 1 X c 1 
get_nodes 
I nver5e_c 
g e t _ n d e s 
node_col 
bl each( xc 
get _ nodes 
upd_prece 
If ( ( desc 

precendants. c 
If ( ( p r e c 

precendants. c 



ched_col ors; 
I nt X, I nt Y) { 
e_c nt + + ; 
r a 1 r It; 
ymap; 

ol or_nodes( ) ; 

or s=[node_col ors . count ( X) >0) ?cn[ X] ; empt ymap; 



h() ; 



rs. erase(X); 
1 r s ) ; 
0; 

n d a n t s ( ) ; 

e nda nt s . count ( St a r t nod e ) ==0) | | ( de s c e n d a nt 5 [ s t a r t node] . si ze( ) ==0) ) {descendant s . cl ear () ; ^ 

learO;} 

e nda nt s . c ount ( t op n od e ) ==0 ) | | ( pr e c e n d a n t s [ t p node 1 . s I z e ( ) ==0) ) {des c e nd a nt s . c I ea r ( ) ; ^ 

learO;} 



void bleach! ma pl& xcolors) 
mapl::lterator It; 
ma pi nO=nodes; 



for ( I t =nO. begi n( ) ; I t ! =nO. end( ) ; I t ++) { 
Int node = (int)lt->flrst; 
mapl::iterator jt; 
dep_col ors=xcol ors; 

get dep col ors(cn[ node] , xcol ors, col or_hi erarchy, 
get _dep_ c ol or s ( c n[ node] , xc ol or s , 1 c h, ly ; // up 
mapdi ff( cnl node] , depcol ors) ; 
remove_from_col or_hi erarchy(dep_col ors); 
map diff(node colors[node],dep colors); 
if [cn[ node] . ii ze( ) ==0) { 



node 



1) ; / / going down 



get _ I 
ma pi: 
for I 



col ors 



nda nt 


s . erase 


( node) 


color 


s . erase 


( node) 


s ( n d 


e , d es c e 


nd a nt s 


s ( n d 


e, prece 


nda nt s 


odes ( 


) ; 




liter 


a 1 r c 1 




c i =c n 


[node]. 


b e g 1 n ( 


nt c = 


( 1 nt ) cl 


■ >f 1 r s 


ol or _ 


nodes . e 


rase( c 


ect ed 


_col ors 


0; 



) ; cl! =cn[ node! . end( ) ; cl ++) { // for node we remove - deactivate I t ' s 



} 

vol 



} 

clear 



d cl ear_di sconnected_col ors( ) { 
If ( des c endant s . s i ze( ) ==0) return; 
I nverse_ch( ) ; 
ma p 1 r ml ; 

map2::iterator it; 

for ( i t =i c h . beg i n ( ) ; i t I =i c h . e n d ( ) ; 



it++) { 
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i nt c = ( i nt ) i t - >f i r s t ; 

if ( ( i chl c] . s i ze( ) ==0) &&( t op_col or 5 . count ( c ) ==0) ) { 
r ml [ c] =t r ue; 

} 

} 

remove_from_col or_hi erarchy(rml ); 
i nverse_ch(I; 

id remove_from_color_tiierarciiy(mapl& x) { 
mapl::iterator it; 

for ( i t =x. begi n( ) ; i t I =x. end( ) ; i t + + ) { 
i nt c = ( i nt ) i t - >f i r St ; 
if (lactive_color(c)) { 

coior_hierarchy.erase(c); 
ma p 1 : : i t e r a 1 r it; 

for ( i t =i cli[ c] . begi n( ) ; i t I =i cti[ c] . end( ) ; i t + + ) { 
int sc = (int)it->first; 
color tilerarciiy[sc].erase(c); 

} 

cont I nue; 



} 



} 

remove_one_col or_from_col orhi erarctiy(c) 



tia 



id remove_one_coior_from_coior_hierarctiy(int c) { 

// given coFor may liave several I sup colors: I.e. red is owned by blue and brown in different 

its before merging together 
mapl:;iterator it; 
ma pi X = coior_hierarcliy[c]; 
col or_hi erarcliy, erase(c); 
coior_nodes.erase(c]; 
if ( i ch[ c] . si ze( ) = = 1) { 

1 1 col or_hi erarchy. erase( i ch[ cl. begi n() - >fi rst); 

I ch. erase(c); 

} 

for ( i t =x. begi n( ) ; i t ! =x. end( ) ; i t + + ) { 
int s c = ( 1 n t ) i t ■ >f i r s t ; 
if (!active_color(sc)) continue; 
I f ( sc==empt ycol or ) cont I nue; 
If ( I ch[ sc] . si ze( ) <=1) { 
I ch. er ase( sc] ; 

removeonecol orfromcoi orhl erarchy(sc) 
nee this owner (c) was removed now we will finally delete sc 
} 

} 



// the other owner may be removed later, but 



id extract_path() { 

path,ciear(); 

booi b=has_path_to_top(top_coi ors. begi n() - >fl rst,emptymap,coi or_hi erarchy, -1); 

coiorpath=path; 

if ( ! b) { 

printf(" No path to top \n"); 

} 

transiate_to_nodes(); 
int pathien=pat[i.size(l; 
if ( pat hi en==gmax_sl ack) { 

printf(" Path Is confirmed \ n"); 
} else { 

prlntf(" Path Is not confirmed \n"); 

} 

id t r a ns I a t e_t o_ nodes ( ) { 
ma pi node_path; 
mapl::lterator It; 

for ( I t =pat h. begi n( ) ; I t ! =pat h. end( ) ; I t ++) { 
Int c= (int)lt->flrst; 
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node_pat h[ col or_nodes[ c] 1 =t rue; 

p r i n f f ( " path %d %d \ n " , c , c 1 r n d e 5 [ c ] ) ; 

} 

path=node path; 

} 

ma pi dep_colors; 

void get_dep_colors(mapl& colors, ma pl& xcolors,map2& fiber, Int dir) { 

// traverse to top for each color in colors avoiding xcolors, those that dont have path to top are 1^ 
d e p colors 

mapl::iterator it; 

for ( 1 t =c 1 r s . be g 1 n ( ) ; 1 t ! =c 1 r s . en d ( ) ; 1 t ++) { 

int c = ( i nt ) 1 t - >f 1 r s t ; 

if ( ( c = = e mpt y c 1 r ) &&( d 1 r ==- 1 ) ) continue; 
path.clear();visited.clear(); 

If (lhas path to topic, xcolors, fiber, dIr)) dep colors[cl=true; 
} ' ' " 

} 

ma pi path, colorpath, visited; 

bool has_path_to_top(int src,mapl& xcolors, ma p2& ch,int dIr) { // from src avoiding xcolors through ch 
mapl::lterator It; 

if ( ( dl r = = 1) &&( chl src] . si ze( ) ==0) ) return true; 
for ( I t =ch[ s r c] . begl n( ) ; I t ! =chl s r c 1 . end( ) ; 1 t ++) { 

int c = (lnt)it->first; 

if ( I a c t i V e_ c 1 r ( c ) ) continue; // Inactive color 
node is removed color becomes inactive 
if ( X c 1 r s . c u nt ( c ) >0 ) continue; 
if ( pat h. si ze( ) >0) { 
If ( dl r = = 1) { 

if (lmap_contains(subcolors[c],path)) continue; 
} el se { 

if ( ! c in every sub( c) ) cont 1 nue; 

} 

} 

path[c]=true; 

if ( ( ( ch[ c ] . si ze( ) = = 0) &&( dl r = = 1) ) I I ( ( c= = empt ycol or ) &&( dl r ==■ 1) ) ) return true; 
bool ret =has_path_to_top(c, xcol ors, ch, dl r); 
if (ret) return true; 



every color is cn[] for so me node, If that 1^ 



} 

path. erase(src) ; 
return false; 

} 

bool c_l n_every_sub( Int c) { 
mapl::lterator It; 

for ( 1 t =pat h. begi n( ) ; 1 t I =pat h. end( ) ; 1 t ++) { 
int sc = (int)it->first; 
if ( subcol or si sc] . count ( c) ==0) { 
return false; 

} 



return true; 



int depth; 

ma pi get_colored_nodes(lnt color) { 
ma p 1 ret; 



ma p 1 
for 



} 



literator it; 

t=nodes.begin();itl=nodes.end();it+ + ) { 
nt node = (lnt)it->first; 

f ( node_col ors[ node] . count ( col or ) >0) r et [ node] =t r ue; 



return ret; 

} 

ma pi traversed; 

I 1 St <l nt > f ul I _pat h; 

ma pi traversed_color_nodes; 

ma pi sup_colors; 

ma p 1 e mp t y ma p ; 
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}; 

struct base_edg : dg { 

static const i nt niax_ nodes =1000; 

bool **docked; // dock registry at the node 

base_edg( ) { 

docked = new booi*[niax_nodes]; 
for ( i nt k=0; k<max_nodes; k+ + ) { 
docked[k]=new bool[max_nodes]; 

for ( i nt j =0; j <max nodes; j ++) d c k e d [ k 1 [ j 1 =f a I s e ; 

} 

} 

virtuai bool add delete(int node, i nt X)=0; 

}; 

struct base_bot : utii { 
int current_node; 
booi done; 
base edg * p; 
base Jot ( ) {}; 

base_bot(int node, base_edg *p_) : c u r r ent _node( node ) , p(p_){}; 
virtuai void advance!) =0; 
virtual bool has_botlets()=0; 

}; 

struct bot : public base_bot { 

map<int, bot> botlets; // newborn bots 
botO {}; 

bot(int node, base_edg ♦p_) : base_bot ( node, p_) { 
done=f al se; 

}; 

voi d advance! ) { 

botlets. clear!); 

ma pi 5id=p - >descendants[ current_nodel ; 
ma p 1 : : 1 t e r a 1 r mi t ; 
int outbounds=0; 
d n e =t rue; 

_ASSERTE! _CrtCheckMemory! ) ); 

int e n t r y _ n d e =c u r r e nt _ nod e ; 

for I mi t =d. begi n! ) ; mi t I =d . e n d I ) ; mi t ++) { 

int node_out =mit->first; 

//if I vi si ted. count ! node_out ) ) continue; 

if I I p- >add_del et e! ent r y_ node, node_out ) ) continue; 

if I I p- >docked[ node out 1 fst epl ) { 
p->dockedInode_out][stepl=true; 

if ! out bounds >0) bot I et s [ bot I et s . s i z e I ) ] =bot | node_out , p) ; 
else { 

current_node=node_out; 

done=f aFse; 

} 

out bounds ++; 

} 

//_ASSERTE! _Crt CheckMemoryl ) ); 

if ( ! n d e_ u t = = e n d _ n d e ) &&I u t bo u n d s = = 1 ) ) done=true; 

} 

if !outbounds= = 0) done =t rue; 

} 

bool has_bot I et s I ) { 

return | bot I et s . s i ze! ) >0) ; 

} 

}; 

struct slacks { 
int this_node; 
bool act 1 ve; 

map<int,dg> *sgs; // dg per slack 
slacks!) {}; 

siacksS operator =!const slacksS s2) { 
this_node=s2.this_node; 



APPENDIX 



18 



C:\mydoc5\deform_\sd9.h 



acti ve=s2. acti ve; 
sgs=new map<i nt , dg>( ) ; 
* s g s = * s 2 . s g s ; 

return *this; 

} 

5lacks(const slacks& s2) { 
thi s_node=s2. thi s_node; 
acti ve=s2. acti ve; 
sgs=new map<int,dg>(); 
*sgs = *s2.sgs; 

} 

siacks(int this_node_) t h i s _ nod e ( t h i s _ n d e_ ) { 
sgs=new niap<int,dg>(); 
a c t i V e =t r u e ; 

}; 

void add_siack(int new_node, int prev_node) { 
thi s_node=new_node; 
if ( sgs- >si ze( ) >0) { 

for (int 5 1 a c k =g ma x_ s I a c k ; s i a c k >0; s I a c k- - 
if ( ( * s gs ) . c ount ( 5 1 ac k) ==0) continue; 
dg sdg = ( *sgs) [ si ack] ; 
if ( sdg . des cenda nt 5 . s i z e( ) ==0) { 
sgs - >er as e( 5 1 ac k) ; 
cont i nue; 

} 



s dg 


. si ack_i evei 


=s 


1 ack+1 




s d g 


. add( new_nod 


e , 


pr ev n 


ode) ; 


( *s 


gs).erase(si 


a c 


k); 




s dg 


. get_nodes(l 








i nt 


n u m_ n d e s =s 


dg 


. nodes 


.size 


( *s 


gs) [ sdg. si ac 


k_ 


i evei 1 


=sdg; 


p r e 


v_ n de ==s t a r 


t_ 


node) 


{ 


dg 


sdg; 








s dg 


. add(new_nod 


e , 


p r e v_ n 


ode) ; 


sdg 


.slack i evei 


= 1 






( 's 


gs)[l]=sdg; 









e i s e 
i f 



} 
} 

} 

siacksS operator +=(slacks& s2) { 
// mt aching dgs are combined 
ma p <i nt , d g >: : i t e r a 1 r mi t ; 

for ( mi t =s g s ■ >b e g i n ( ) ; mi t ! =s g s - >e n d ( ) ; mi t + + ) { 
int slack_level=(* mi t). first; 
if ( s2. sgs- >count ( si ack_l evei ) >0) { 

// merge the matching ones: 

(*sgs)[slack_level]+=(*s2.sgs)[slack_ievel 



■>end(); mi t + + ) { 

[ si ack i evei ] ; 



// 


new 


ones are added 




for 


( mi 


t=s2.sgs->begin(); 


mi t I =s 2 . s g s 




i nt 


si ack_l evei =(*mi t) 


.first; 




i f 


(sgs->count(siack_i 


evei)= = 0) { 






( *sgs) [ si ack_i evei 


l=(*s2.sgs) 



} 

return *this; 



booi rm_node(int X,int Y) { 
ma p <i n t , d g >: : i t e r a t r mi t ; 
i i St <i nt > r ml ; 

for ( mi t =s g s - >beg I n ( ) ; mi t ! =s g s - >e n d ( ) ; mi t + + ) { 
dg &g=( * mi t ) . second; 
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i n t slack = ( * mi t ) . f i r s t ; 

if (( g. pr ecenda nt s . count ( X) ==0) &&( g. des c endant s . c ount ( X) >0) ) { 
pr i nt f ( " WTF" ) ; 

} 

if ( g. pr ecendant 5 . count ( X) ==0) continue; // one may bypass 4 otiiers go through 4, stiii need to^ 
process them !!! 

g.rm_node(X,Y); 

g.clear_inacti«e_colors(); 

g.slack_level=(* mi t). first; 

( * sgs ) [ ( * mi t ) . f i r St ] =g; 

if (g. descendants, si ze() ==0 ) { 

rmi.push_baci(((*mit).first); 

c n t i n u e ; 

} 

} 

iist<int>::iterator it; 

for ( i t =r mi . be g i n ( ) ; i t I =r mi . e n d ( ) ; i t + + ) { 

sgs->erase((int)(*it));// iose that dag compieteiy 

} 

_ASSEf!TE( Crt Checl(Memory( ) ); 
return ( r mi . s i z e ( ) >0 ) ; 

} 

}; 

struct edg : base_edg { // extended directed graph - slacl< graph is attached per node and siack ievel 
siacks *node_slacks; // int node mapped onto a slacic map which stores sg for each siack 
edg( ) {hpat h^cnt =2; } 
int siack_depth_c[ieck; 
int e i ; 

map<int,bot> bots; 
int hpath_cnt; 

booi a d d_ d e i e t e ( i n t node, int X) { 
add_deiete_cnt+ + ; 

siacks ns(node_slacks(node]); // make a separate copy for each descendant 

if (ns.rm_node[X,node)) { // if this siacks stack sgs which contain node "it", they need to be ^ 

pruned 

if ( ns. sgs- >si ze( ) ==0) return faise; 

} 

ns.add_slack(X,node); 

if ( node_sl acks[ X] . act i ve) node_s i ac ks [ X] +=ns ; 
ei se node_sl acks[ X] =ns; 
return true; 

} 

void bot_march() { // makes only one step at a time, returns control back 
ma p 1 visited ; 

hot s [ 0] =bot ( St ar t _node, this); 
keyl=0; 

booi done=faise; 

for(step = l;bots.size()>0;step+ + ) { 
ma p <i nt , bot >: : i t er at or it; 

ma p 1 r mc ; 

map<int,bot> new_bots; 
if ( St ep>gmax_sl ack) { 

printf(" step>gmax siack \n"); 

} 

for ( i t =bot s . beg i n ( ) ; i t I =bot s . en d ( ) ; i t ++) { 
int k e y =i t ■ >f i r s t ; 
i t - >second. advance! ) ; 

if ( node_sl acks[ end_node] . sgs- >count ( gmax_sl ack) >0) { 
bot s . c I ea r ( ) ; 

p r i n t f ( " ma r c h c mp I e t e d \ n " ) ; 

done=t r ue; 

break; 

} 

if ( I i t - >second. done) { 

i f (i t->second. has_boti ets() ) add_maps<map<i nt, bot>>( &new_bots, &i t - >second. boti ets); 
} else rmclkey]=true; 
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} 



} 

i f ( done) break; 
ma p 1 : : i t e r a t r mi t ; 

for ( mi t =r mc. be g i n ( ) ; mi t ! =r mc . e n d ( ) ; mi t + + ) b ot s . e r a s e ( mi t ■ >f i r 5 1 ) 
add_maps<map<int,bot>>(&bots,&new_bots); 

p r i n t f ( " ma r c h c mma n d %d , n u mb e r b 1 s %d \ n " , s t e p , b 1 s . s i z e ( ) ) ; 



} 

void init_siacks() { 

node_si acks=new si acl(s[gmax_si ack+2]; 

n e w_ c i r =1 ; 

ma p 1 : : i t e r a 1 r mi t 2 ; 

ma pi &nodes = get_nodes(); 

//_ASSERTE( _CrtCheckl^emory( ) ); 

for ( mi 1 2 =n d e s . b e g i n ( ) ; mi 1 2 ! =n od es . e n d ( ) ; mit2++) { 
i nt X=( * mi 1 2) . f i r St ; 
node si acksl X] =sl acks( X) ; 

} 



} 

voi 



} 

voi 



d u pd 


St 


art endO { 




// St 


ar 


t node i 5 


ne that 


ma p2 : 


: i 


t er at or it; 




for ( 


i t 


=descendant 


s . beg i n 


i 


nt 


node = i t ■ 


>f i r s t ; 


i 


f 


( precendant 


s. count 






St art node 


=node; 






break; 




} 








} 








for ( 


i t 


=precendant 


s . b e g i n 


i 


nt 


node = i t - 


>f i r s t ; 


1 


f 


( descendant 


s . count 






end_node=n 


ode; 






break; 




} 








} 








d r u n i 


t( 


i nt 1) { 




5 u be 


i ors. cl ear ( ) ; 




ei =1 ; 









has no pr ecendant s 
( ) ; 1 t ! =descendant s . end( ) ; i t ++) { 
(node)==0) { 

(); 1 t ! =pr ec endant s . end ( ) ; i t ++) { 
(node)==0) { 



//_ASSERTE( _CrtCheckMemory( ) ); 
i oad gr aph( 1 ) ; 
upd_start_end( ) ; 

init_siacks(); 
b 1 _ ma r c h ( ) ; 

if ( ( * n d e_ s I a c k s [ e n d _ n d e ] . s g s ) . c u n t ( g ma x_ s I a c k ) >0 ) { 
pr i nt f ( " H- Gr a ph %d \ n" , i ) ; 

(*node_slacks[end_node],sgs)[g ma x_siack], extract path(i 
} ei se { " 

pr i nt f ( " Not an H- Gr aph %d \ n" , i ) ; 

} 
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#i 


nc 


ude 


<i ost r eam> 


#i 


nc 


ude 


<f St r eam> 


#i 


nc 


ude 


<l i St > 


#i 


nc 


ude 


<s t r i n g > 


#i 


nc 


ude 


<d e q u e > 


#i 


nc 


ude 


<l i s t > 


#i 


nc 


ude 


<map> 


#i 


nc 


ude 


<t i me . h > 



using na me space std; 

#def i ne smal I num - 99 9 9 9 9 9 

#def i ne I argenum 99 9 9 9 9 9 

typedef map<int,bool> ma pi; 
typedef map<int,mapl> map2; 
int gmax_slack; 



#def 1 


ne 1 s n 


#def i 


ne s 1 g 


#def i 


ne r ou 


St at i 


c cons 


St r uc 


t ut i i 


virtu 


a i voi 


voi d 


LoadCo 


clia 


r * i i n 


i f s 


t r e a m 


f 5 . 


p e n ( " 



n(x) ((X) != (X)) 

(X) (((x)>0)?l:(((x)= = 0)?0:-l)) 

d( X) ( ( ( X) - f I oorl xl >cei I ( X) - ( X) ) ?cei I ( X) : f I oor( x) ) 



{ 

d pr oces 
nfi g() { 
e =( c h a r * 



f ( f s. good( ) ==0) 
p r i n t f ( " \ n Error 



} 

char* tag=(c[iar*) 
char* vaiue=(char 
for (int r w = ; f s 

fs.getiine(iine, 

s t r c a t ( I i n e , " \ " 

if (strstr(iine, 

strcpy(tag, strtok( 

strcpy(val ue, strtokf 

processTagVai ue(tag,val ue); 

//// ASSERTEI Cr t CheckMemory( ) ) ; 

} 

f s. ci ose( ) ; 



r aphpat h[ ] = 


" c : W g r a p h s \\ 




sTagVal ue( cha 


r * t a g , c h a r * 


val ue) {} 


) ma 1 1 c ( 1 2 4 ) 






X t " , i s : : i n ) ; 






{ 






r opening the 


c nf i g file. 


ex i t i n g . \ n 


opening the 


c nf i g f i i e , 


e X i t i n g . \ n " 


ma 1 1 c ( 2 5 6 ) ; 






* ) ma 1 1 c ( 2 5 6 ) 






.good()l=0;ro 


W + + ) { 




, 1 2 4, ' \ n' ) 






"); 

, " =" ) = = NULL) 


c nt i n ue ; 




ok( i i ne, " =" 


)); 




rtok( NULL, " 


=" )); 





static const 


i nt 


ms z 


=1 


int s z i ; 








map2 descenda 


nts; 


// 


n 


map2 precenda 


nt s ; 


// 


f 


booi ioad gra 


ph( i 


nt 


i ) 


/ / M = [ 


1 


1; 





to somewh 


ere, 


i . 


e . 


/ / load d 


e s c e 


n d a 


nt 


descendan 


t t h 


ose 


a 


char line 


[ ms z 


1 ; 




szi =( 1 nt ) 


r oun 


dis 




i f s t r ea m 


f s; 







char f n a me [ 6 1 ; 

s p r i n t f ( f n a me , " %s g r a p h %d . t X t " , g r a p h p a t h , i 

fs.open(fname); 

i f (f s. goodO ==0) { 

f s . c i OS e( ) ; 

return false; 
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} 

fs.getline(line, msz, '\n'); 

5 1 r c a t ( I i n e , " \ " ) ; 

char * pcti = strtok (line, ";[]"); 

int nn = msz/szl, k=0; 

char **node=new char*[nn+l]; 

for ( k=l; pch ! = NULL; k+ + ) { 

p c h = s t r t k (NULL, ";[]"); 
f (pch= = NULL) break; 

node[k]=new char[szl]; 

strcpy(node[k],pch); 

s t r c a t ( n d e [ k ] , " \ " ) ; 

} 

for (int j =1; j <k; j + + ) { 

char * pch = s t r t o k ( n o d e [ j ] , " " ) ; 
for ( i nt ni =1; pch ! = NULL; ni + + ) { 
int val = atoi(pch); 
if (val = = 1) { 

de5cendantslj][ni]=true; 
if (j ==ni ) { 

pr i nt f ( " j ==ni \ n" ) ; 

} 

i f ( descendant s . count ( nl ) >0) { 

if ( des cenda nt s [ n i ] . c u nt ( j ) >0 ) { 

p r i n t f ( " circular reference \ n " ) ; 

} 

} 

precendantslni][j]=true; 

} 

pch = St rt ok (NULL, " " ) ; 
If (pch= = NULL) break; 

} 

} 

gmax sl ack=de5cendants. si ze( ) ; 
f 5 , c I 05 e( ) ; 
return true; 

} 

ma pi flatten(map2& a) { 
ma p 1 ret; 

map2::lterator it; 

for ( i t =a. begi n( ) ; i t ! =a. end( ) ; i t + + ) { 
Int key = ( i nt ) I t - >f i r s t ; 
r et [ key 1 =t r ue; 
mapl::lterator jt; 

for ( j t =a[ key] . begi n( ) ; j t I =a[ key] . end( ) ; j t ++) { 
int V = (int)jt->flrst; 
r et [ V ] =t r ue; 

} 

} 

return ret; 

} 

void gen_graph(lnt I) { 

Int num_nodes=17,locai_degree=2; 
while (true) { 
map 2 rows; 

map<lnt,lnt> node_lnds; 
ma p 1 c h k ; 

bool restart =faise; 

for (int k=l; k<num nodes-1; k+ + ) {111 and num nodes are forbidden 
f or ( i nt 5=0; t r ui; 5 ++) { 

int j =r a nd ( ) %n u m_ nodes ; 
if ( 5>5*num_nodes) { 

restart =t rue; // one last available number may match the index variabie, restart then 
break; 

} 

if ( ( chk. count ( j ) >0) I I ( j <2) I I ( j = = k) ) continue; 
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node_inds[kl=j; 
c hk[ j ] =t r ue; 
break; 

} 

if (restart) break; 

} 

if (restart) continue; 
map<int,int> path; 
int node_ind=0; 
i nt pr ev_ r =1 ; 

for (int r=l; r <n u m_ n o d e s ; r + + ) { // n-1 rows 
if ( n d e_ i n d s [ r ] = = ) nod e_ i n d s [ r ] =nu m_ n o d e s ; 
if (prev_r==node indslr]) { 
i nt i =0; 

} 

r ws [ pr e v_ r ] [ node_ i nds 1 r ] ] =t r ue ; // node r connects to node_ind 
if (path.caunt(prev_r)>0) { 
i nt i =0; 

} 

patti[prev_r]=node_inds[r]; 
pr ev r =nQde indslr]; 

} 

for (int r=l; r <n u m_ nodes ; r + + ) { // n-1 rows 

if ( ( r >1) &&( r <num_nodes- 1) ) { // one before last only has last one as descendant; 
i nt j =1; 

wh i i e (true) { 

doubl e d=rand( ) / ( ( doubi e) RAND_MAX) ; 
int k=(int)(d*num_ nodes); 

if (k<2) continue; // 1 is start node, it's not anybodies descendant 

if ( k = = n u m_ no d e s ) continue; // num_nodes is end node, it's not anybodies descendant 

if ( r ws [ r ] . c u n t ( k ) >0 ) continue; 

if ( r ws [ k ] . c u n t ( r ) >0 ) continue; // at 6 we have 15 then at 15 we cant have 6 
if ( r = = k ) continue; 
r ows [ r ] [ k] =t r ue; 

j + + ; 

if ( j >i ocal degree) break; 

} 

} 

} 

save_graph(rQws,i); 
char fnameI1281; 

s p r i n t f ( f n a me , " %s h p a t h %d . t X t " , g r a p h p a t h , i ) ; 

save_map( path, fname) ; 

p r i n t f ( " path length %d " , p a t h . s i z e ( ) ) ; 

break; 

} 

} 

void s a V e_ g r a p h ( ma p 2 &rows,int i){ 
char fname[128]; 
int n u m_ n d e s =r ws . s i z e ( ) +1 ; 
s p r i n t f ( f n a me , " %s g r a p h %d . t X t " , g r a p h p a t h , i ) ; 
map2::iterator it; 
string s t r ( " M = [ " ) ; 
int r =1; 

for ( i t =r ws . be g i n ( ) ; i t I =r o ws . e n d ( ) ; i t ++) { 
ma p 1 &r w = i t ■ >s e c n d ; 
for (int k=l; k < = n u m_ nodes ; k+ + ) { 
if (k>l) str+ = " "; 
if ( r ow. count ( k) >0) { 

St r + = "1" ; 
} else St r +=" 0" ; 

} 

r ++; 

if ( r ! =r ws . s i z e ( ) +1 ) s t r +=" ; " ; 

} 

St r + = " 1 " ; 
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ofstream ofs(fname, 
of s << s t r << " \ n " ; 
ofs. close!); 



out 1 



} 

voi 



d 5 a V e_ ma p ( ma p <i nt , i nt > Sia,char *fname){ 
ofstream ofs(fname,ios::out); 
ma p <i n t , i nt >: : i t e r a 1 r it; 
for (it=a.begin(); it!=a.end(); 
of s << i t ■ >f i r St << " - > " << 

} 

of s . c I OS e( ) ; 



it++) { 
i t - >second 



nt get_cols(char *name) { 
char I i ne[ 1 2 4 ] ; 
i f s t r e a m f s ; 
f s . p e n ( n a me ) ; 
if (f s. good!) ==0) { 

f s . c I OS e( ) ; 

return ■ 1 ; 

} 

cliar *token;int col=0; 
f s . get i i ne( 11 ne, 1 24, ' \ 
St r cat ( i i ne, " \ 0" ) ; 
token = strtokl iine, "," 



for(coi=0; token ! 
token = strtokl 



= NULL; 
NULL, 



coi ++) { 



} 



f s . c i s e ( ) ; 
return coi; 

nt get_iines(const char ♦name) { 
char Fine! 1 2 4 ] ; 
i f s t r e a m f s ; 
fs. open( name) ; 
if (fs. good( ) ==0) { 

f s . c i OS e ( ) ; 

return ■ 1 ; 

} 

i nt r ow=0; 

for ( r ow=0; f s . good( ) ! =0 
f s . get i i ne( I i ne, 1024 
if ( St r i en( 

} 

f s . c i OS e( ) ; 
return row; 



row++) { 
'\n' ); 
line)<l) return row; 



} 



tempiate<typename T> 
T* join_maps(T *a, T * 
T : : i t e r a 1 r mi t ; 
for ( mi t =b- >be g i n ( 
i n t k e y =( * mi t ) 



b) { 



} 



i f 



} 



mi t I 

r St ; 

>count(key)= = 0) 



=b->end(); mi t ++) { 



{ 



(*a)[key]=(*b)[key]; 



return a ; 

} 

1 nt keyl; 

template<typename T> 
T* add.mapsIT 'a, T *b) { 
T : : 1 t e r a t r mi t ; 
1 nt j =a- >si ze( ) ; 
for ( mi t =b- >be g 1 n ( ) ; mi t ! : 
i nt key =( • mi t ) . f i r s t ; 
(•a)[keyl]=(*b)[key]; 



b->end(); mi t ++) { 
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keyl++; 

} 

return a ; 

} 

void map difflmapl &a,mapl &b) { // need to avoid bleaching upstream from precendant of X, otherwise we 
wi i i aiways bieach all descendants of X 
mapl::iterator it; 

for ( i t =b. begi n( ) ; i t ! =b. end( ) ; i t + + ) { 
i n t k e y =( * i t ) . f i r s t ; 
if ( a. count ( key) >0) a. erase! key) ; 

} 

} 

booi map_contains(mapl &a,mapl 5ib) { // ma pi has map2 
ma p 1 : : i t e r a 1 r it; 

for ( i t =b. begi n() ; i t ! =b. end( ) ; i t + + ) { 
i nt key =( * i t ) . f i r s t ; 
if ( a. count ( key) ==0) return false; 

} 

return true; 

} 

void merge_maps(map2& a, map2& b) { 
ma p 2 : : i t e r a 1 r mi t ; 

for ( ml t =a . b e g i n ( ) ; mi t ! =a . e n d ( ) ; mi t ++) { 
int node =(* mi t). first; 
if ( b. count ( node) >0) { 

join ma ps(&a [node], &b [node]); 

} 

} 

// those that were missing in this graph just get added fromg2: 

map2::iterator it; 

for ( i t =b. begi n( ) ; i t I =b. end( ) ; I t ++) { 
int node =(int)(*it). first; 
If (a.count(node)= = 0) { 

if ( b. count ( node) >0) a[ node] =b[ node] ; 

} 

} 

} 

void me r g e_ ma p s ( ma p <i nt , i I s t <l nt >>& a, map<i nt , i i st <l nt >>& b) { 
map<int,llst<lnt>>;:lterator mit; 

// those that were missing In this graph just get added fromg2; 
for ( ml t =b. begi n( ) ; ml t ! =b . e n d ( ) ; mi t ++) { 

Int node =( I nt ) ( * ml t ) . f I r s t ; 

If ( a. count ( node) ==0) { 

If ( b. count ( node) >0) a[ node] =b[ node] ; 

} 

} 

} 

void me r g e_ ma p s ( ma p <i n t , i n t >& a, ma p <i n t , i n t >& b) {// color_nodes 
map<int,int>::iterator it; 
for ( i t =b. begi n() ; i t I =b. end ( ) ; I t + + ) { 
int c i r =( I nt ) i t - >f i r s t ; 
a[colorl=b(color]; 

} 

} 

booi map_equais(map2& a, map2& b) { 
return a in b(a,b) &&a In b( b, a ) ; 

} 

bool a_ln_b(map2& a, map2& b) { 
map2::lterator It; 

for ( I t =a. begi n( ) ; I t I =a. end( ) ; I t + + ) { 
int key = (lnt)it->first; 
if ( b. count ( key) ==0) return false; 
mapl;;lterator ]t; 

for ( j t =a[ key] . begi n( ); j t ! =a[ key] . end( ); j t ++) { 
int key2 = (int)jt->flrst; 
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if ( b[ key] . count ( key2) ==0) return false; 

} 

for ( j t =b[ key] . begi n( ) ; j t ! =b[ key] . end( ) ; j t ++) { 
int key2 = (int)jt->fir5t; 
if ( a[ key] . count ( key2) ==0) return false; 

} 

} 

return true; 

} 

ma pi 1 nt e r s ec t ( ma pi &a, ma pi &b) { 
// those in a that In b 
ma p 1 ret; 

ma p 1 : : i t e r a 1 r it; 

for ( i t =a . beg i n( ) ; i t I =a . end ( ) ; i t + + ) { 
int key = (int)it->first; 
if ( b. count ( key) >0) { 
r et [ key] =t r ue; 

} 

} 

return ret; 

} 

}; // uti i 
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#i nc I ude " s t daf x . h" 

#i nc I ude " s d g . h " 

voi d gen_graphs( i nt n) { 

e d g g r ; 

s r and ( t i me( NULL) ) ; 

for ( i n t g =1 ; g <=n ; g ++) { 

p r i n t f ( " generating graph %d \ n " , g ) ; 

gr.gen graph(g); 

} 

} 

i nt _tmain(int argc, TCHAR* argv[l) { 
if (fal se) { 

gen g r aphs ( 1 ) ; 
} else { 

int gi=_ttoi(argv[l]); 

p r i n t f ( " Processing Graph %d \ n " , g i ) ; 

edg gr ; 

gr.runit(gi); 

} 

return ; 



